arXiv: 1503.02061vl [physics.acc-ph] 6 Mar 2015 


A Comparison of Storage Ring Modeling with COSY INFINITY, ZGOUBI, and 
MADS 


Robert Hippie and Martin Berz 

Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA 

Currently there is significant interest in the use of storage rings to search for an electric dipole moment 
(EDM) in hadrons [1]. This requires utilizing the storage ring as a precision measuring device [2]. Part of 
understanding the detailed behavior of storage rings comes from careful analysis of fringe fields [3], but 
the various tracking codes available differ in their ability to model such behavior. It is the purpose of this 
paper to investigate these differences. 

A major storage ring facility actively engaged in the search for hadron EDMs is the COoler SYnchrotron 
(COSY) at Eorschungszentrum liilich [4]. We modeled a simplified version of this storage ring using 
three well-known simulation codes - MADS [5], ZGOUBI [6] and COSY INEINITY [7]. MADS is a 
“transfer map” code of order 2, which means that the state of the particle in phase space is maintained as 
a vector, and the differential equations governing the motion of particle through the storage ring elements 
are represented by transfer maps. To track a particle through a system, one merely needs to perform map 
composition. MADS also has the capability to track particles symplectically using generating functions of 
third order [S]. 

ZGOUBI does not use the transfer map technique, but rather integrates the Lorentz equation by time step¬ 
ping based on a Taylor series in path length. The coefficients of the Taylor series in time are determined by 
an additional Taylor expansion of the magnetic field, to fifth order maximum, if the fields are given analyt¬ 
ically, and by an out of plane expansion based on numerical differentiation otherwise. ZGOUBI has few 
programming capabilities beyond a simple looping mechanism to provide multiple passes through an opti¬ 
cal system. ZGOUBI provides support for fringe fields via Enge coefficients [9]. The software distribution 
also includes a powerful post processing module called ZPOP for plotting and data visualization. 

COSY INFINITY is a combination of the advantages of the transfer map approach and integration codes. It 
is primarily a transfer map code, but utilizes integration internally to create highly accurate maps for fringe 
fields [10]. Built into COSY INEINITY is an interpreter for the specialized COSYScript programming 
language [11], which allows the researcher to simulate charged particle optics systems to a high degree 
of accuracy using the techniques of Differential Algebra. Eringe fields are specified by Enge coefficients 
which can be input by the user to model actual field measurements, or a default set of typical values can 
be chosen. 

To establish a baseline, we begin with a simplified hard-edge model of the COSY storage ring (Eigure 
[^. The ring is highly symmetric, incorporating two 40 m telescope regions, and two 52 m arcs. Each 
arc is composed of three identical bending segments, each with mirror symmetry. The bending elements 
are rectangular dipoles. There are 16 sextupoles (not indicated on the figure) at various locations around 
the lattice. Our model incorporates only the bending and focusing elements - the sextupoles in the actual 
lattice are not modeled. After implementing the storage ring elements into the three codes, we confirm 
that the first order transfer matrices are essentially identical (Table [^. 
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Table I. First order transfer matriees for the three eodes without fringe fields. Note that COSY INFINITY 
and MADS use transfer matriees and thus naturally agree to high aeeuraey, whereas ZGOUBI ealeulates 
the transfer map as a result of integration of nearby orbits whieh is slightly less aeeurate. 

Having verified that our lattiees agree to first order, we perform some initial traeking runs and eompare the 
results. Figureshows the traeking pietures for a single 970 MeV/e proton over 2000 turns with eosine- 
like initial eonditions at transverse amplitude < 1 em. There is quite good agreement between MADS 
and COSY INFINITY when performing traeking at MADS’s maximum order of 2, ineluding the inward 
spiral eharaeteristie of non-sympleetie traeking. The output for ZGOUBI also shows a small violation of 
sympleetieity, but apparently ZGOUBFs aeeuraey is higher than that of seeond order transfer maps. 

With COSY INFINITY and MADS we ean enable sympleetie traeking [12], whieh ZGOUBI does not 




Figure I. Simplified model of the COSY storage ring at Forsehungszentrum liilieh. 










Figure 2. Second order tracking in the horizontal plane created by COSY INFINITY (left) and MADS 
(center), compared with ZGOUBI (right) which does not use transfer maps, under identical initial con¬ 
ditions and lattice parameters. No symplectification is used, showing that second order is insufficient to 
describe the dynamics. 



Figure 3. Varying transverse initial conditions from 1 cm to 9 cm in the horizontal plane for second 
order tracking with COSY INFINITY (left) and MADS (center) and compared with ZGOUBI (right). 
Symplectic tracking is enabled in COSY INFINITY and MADS. This is not available in ZGOUBI, which 
shows a widening of the orbit bands indicative of violation of symplecticity. 



Figure 4. Same tracking as Figure now utilizing COSY INFINITY’S 9th order tracking without sym¬ 
plectification, still without fringe fields. The high order of the map and resulting accuracy avoids the 
violation of symplecticity visible in Figure 2. Even to this high order, there is very little nonlinearity 
evident in the dynamics. 

support. Figure [^exhibits tracking pictures with symplectic tracking enabled. They are virtually identical 
for COSY INFINITY and MADS, the phase space orbits having a well-defined elliptical shape. This es¬ 
tablishes a baseline between COSY INFINITY, ZGOUBI and MADS. Turning to the question of dynamic 
aperture, we track particle orbits in 1 cm steps from the reference orbit out to 9 cm - twice the physical 
aperture of the actual ring. We see little significant deviation from linear behavior across all three codes. 
So far, with COSY INFINITY and MADS, we are still tracking with second order transfer maps. With 
COSY INFINITY we can push to higher orders with very little extra processing time. Figure|^is the same 
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Table 2. First-order transfer matrices for COSY INFINITY and ZGOUBI with full fringe fields enabled. 
The horizontal trace (the sum of the bold matrix elements) is greater than 2 in magnitude, which indicates 
instability. 

tracking run, this time to 9th order in the transfer maps. Still there is very little nonlinearity evident in the 
dynamics, and the dynamic aperture appears unlimited. 

Judging from these results alone, one might conclude that the lattice is remarkably stable. We now inves¬ 
tigate the effect of incorporating fringe fields. COSY INFINITY has the convenience of enabling a default 
set of fringe field profiles with the single command “FR 3”. Table shows the COSY INFINITY first 
order transfer map with fringe fields turned on. Notice that the motion in the horizontal direction is now 
unstable (ITracel > 2). 

We can confirm this by adding identical fringe fields to ZGOUBI and checking the computed transfer 
matrix (Table |^. A similar horizontal instability is indicated in the ZGOUBI transfer matrix as well. It is 
apparent that some adjustment of the magnet strengths is necessary to stabilize the beam. Without making 
any specific assumptions regarding which adjustments to make, we choose an arbitrary, minimally inva¬ 
sive approach which respects the existing symmetries of the lattice. The numbered elements in Figure 
indicate quadrupoles which are equal in field strength. By preserving this equality, we maintain the sym¬ 
metry that was designed into the lattice. This allows 8 degrees of freedom. Since we are trying to match 8 
elements of the transfer matrix, we require only six degrees of freedom due to the unity determinant. By 
utilizing an implementation of the MINPACK LMDIF optimizer built into COSY INFINITY [13], we find 
a (not necessarily unique) set of quadrupole strength multipliers which achieve our objective: 

Xi = 1.022183, X 2= 1.019331, X 3= 1.009453, X 4=1.019948, 

X5= 0.999410, X6= 1.019427, X7= 1.031603, X8=1.013424. 

We see this requires an adjustment of only a few percent in the quadrupole strengths. Implementing these 
changes into all three codes, we successfully restore the lattice to its original design transfer matrix. 

Figure]^ shows the resulting tracking pictures for COSY INFINITY. This run was 9th order in the transfer 
matrices, with the fringe fields enabled for all lattice elements. In the horizontal direction, nonlinearities 
are readily apparent. The shape of the phase space ellipses are seen to vary with increasing distance from 
the reference orbit. At 8 cm, seven “islands” of stability appear. In the vertical direction, there is substantial 
beam loss at distances > 6 cm. Neither of these phenomena are visible in the hard-edge approximation. 

Figure shows the same tracking experiment using ZGOUBI. It is interesting to note that the seven 






Figure 5. COSY INFINITY 9th order tracking with full fringe field simulation capabilities in vertical 
(left) and horizontal (right) planes. 



Figure 6. ZGOUBI tracking with full fringe field simulation capabilities in vertical (left) and horizontal 
(right) planes. The fringe field Enge coefficients are identical to those used in COSY INFINITY. The 
rough structure and stability boundaries are similar to those of Figure 5, but symplecticity violations are 
apparent. 

“islands” revealed by COSY INFINITY are visible with ZGOUBI as well but these have the appearance 
of “gaps” in the tracking picture. In the vertical direction, ZGOUBI also agrees with COSY INFINITY’S 
prediction of substantial beam loss at distances greater than 6 cm. 

In conclusion, modeling a storage ring without taking fringe fields into account provides an overly op¬ 
timistic dynamic aperture and an incomplete picture of the dynamics. Although COSY INFINITY and 
ZGOUBI differ in their methods of tracking, both codes hint at similar nonlinear dynamics and agree very 
well for the linear motion with fringe fields. This agreement across codes, despite completely different 
tracking methods, gives confidence that the predictions made reflect real physics. There is a substantial 
difference in run time, however (Table|^, and for moderately large numbers of turns and particles, COSY 
INFINITY is about three orders of magnitude faster than the identical simulation in ZGOUBI. MADS, 
which does a good job of modeling within the SCOFF (Sharp Cutoff Fringing Field) approximation, has 
limited fringe field capabilities. This limits its usefulness for lower momentum tracking over many turns, 
precisely the type of tracking required for mid-size storage rings [14]. 
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Table 3. Comparison of tracking execution times of COSY INFINITY and ZGOUBI at their respective 
maximum preeisions (order 9 transfer map for COSY INFINITY, 5th order Taylor series integration for 
ZGOUBI). ZGOUBI exeeution times are proportional to the number of turns and around 0.187 seeonds 
per turn per partiele. COSY INFINITY requires an initial investment in the eomputation of a transfer map, 
but for larger turn numbers tracks for 1/50000 of a second per turn. 
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